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Abstract 

A novel method is proposed to determine an analytical expression for eigenfunctions and numer- 
ical result for eigenvalues of the Schrodinger type equations, within the context of Taylor expansion 
of a function. Optimal truncation of the Taylor series gives a best possible analytical expression 
for eigenfunctions and numerical result for eigenvalues. 
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I. INTRODUCTION 



One of the source of progress of the sciences depends on the study of the same problem 
from different point of view. Besides their progress of the sciences those different point 
of views include a lot of mathematical tastes. Over the years considerable attention has 
been paid to the solution of Schrodinger equation. Determination of eigenvalues of the 
Schrodinger equation via asymptotic iteration method (AIM) has recently attracted some 
interest, arising from the development of fast computers jl-4]. This method have been 



widely applied to establish eigenvalues of the Schrodinger type equations ilrQ- Although 
the AIM formalism is very efficient to obtain eigenvalues of the Schrodinger equation, it 
requires tedious calculations in order to determine wave function if the system is not exactly 
solvable. When Schrodinger equation includes a non solvable potential, the calculation of 
wave function involved with a large number of terms will lose its simplicity and accuracy. 

In this paper we will discuss new formalism based on the Taylor series expansion method, 
namely Asymptotic Taylor Expansion Method (ATEM). Although, Taylor Series Method 



ll| is an old one, it appears, however, not has been fully exploited in the analysis of both 



in solution of physical and mathematical problems. Yet, even today, new contributions to 
this problem are being made 12|. Apart from its formal relation to AIM, ATEM has also 
been easily applied to solve second order linear differential equations by introducing a simple 
computer program. We would like to mention here that the ATEM is a field of tremendous 
scope and has an almost unlimited opportunity, for its applications in the solution of the 
Schrodinger type equations. One can display a number of fruitful applications of the ATEM 
in different fields of the physics. For instance, our formalism of ATEM gives a new approach 
to the series solution of the differential equations as well as interrelation between series 
solution of differential equations and AIM. The method can also be applied to solve Dirac 
equation and Klein Gordon equation. In this paper we address ourselves to the solution of 
the eigenvalue problems by using the ATEM. 

One of the fundamental advantage of ATEM is that (approximate) analytical expressions 
for the wave function of the associated Hamiltonian can easily be obtained. We note that 
ATEM also gives an accurate result for the eigenvalues when it is compared to AIM. 

It should also be noted here that the determination of wave function by using AIM require 
tedious calculations if the system is not exactly solvable. For non exactly solvable potentials, 
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the calculation of wave function by using AIM involved with a large number of terms will lose 
its simplicity and accuracy. Therefore, the method introduced here is useful to determine 
an analytical expression for the wave function of the non exactly solvable equations. 

The paper is organized as follows. The first main result of the paper is given in section by 
reformulating the well known Taylor series expansion of a function. Section 3 is devoted to 
the application of the main result for solving the Schrodinger equation including various po- 
tentials. As a practical example, we illustrate solution of the Schrodinger equation including 
anharmonic oscillator potential and the Hamiltonian of an interacting electron in a quantum 
dot. In this section we present an approximate analytical expression for eigenfunction and 
numerical results for eigenvalues of the anharmonic oscillator potential and the Hamiltonian 
of an interacting electron in a quantum dot. We also analyze asymptotic behavior of the 
Hamiltonian. Finally we comment on the validity of our method and remark on its possible 
use in different fields of the physics in section 4. 



II. A NEW FORMALISM OF THE TAYLOR EXPANSION METHOD 

In this section, we show the solution of the Schrodinger equation for a quite ample class 
of potentials, by modifying Taylor series expansion by means of a finite sequence instead 
of an infinite sequence and its termination possessing the property of quantum mechanical 
wave function. In quantum mechanics bound state energy of the atom is quantized and 
eigenvalues are discrete and for each eigenvalues there exist one or more an eigenfunctions. 
When we are dealing with the solution of the Schrodinger equation we are mainly interested 
in the discrete eigenvalues of the problem. The first main result of this conclusion gives 
necessary and sufficient conditions for the termination of the Taylor series expansion of the 
wave function. 

Let us consider Taylor series expansion of a function f(x) about the point a: 
fix) = f(a) + (x - a)f(a) + i(x - + \(x - a) 3 f 3 Ha) + • • ■ 

n=0 

where f^ n \a) is the n th derivative of the function at a. Taylor series specifies the value of 
a function at one point, x, in terms of the value of the function and its derivatives at a 
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reference point a. Expansion of the function f(x) about the origin (a = 0), is known as 
Maclaurin's series and it is given by, 

/(*) = /(0) + xf'(0) + ]-x 2 f"(0) + i* 3 / (3) (0) + ■ ■ ■ 

2 o 

oo n 

= E^/ (n) (°)- ( 2 ) 

n=0 

Here we develop a method to solve a second order linear differential equation of the form: 

f"(x)=p (x)f(x)+q (x)f(x). (3) 

It is obvious that the higher order derivatives of the f(x) can be obtained in terms of the 
f(x) and f'(x) by differentiating ([3]). Then, higher order derivatives of f(x) are given by 

f n+2 \x) = p n (x)f(x) + q n (x)f(x) (4) 

where 

p n (x) = po(x)p n -x(x) +p' n _ 1 (x) + q n -x(x), and 

q n (x) = q (x)p n -i(x) +q' n _ l (x). (5) 

Of course, the last result shows there exist a formal relation between AIM and ATEM. To this 
end, we conclude that the recurrence relations ([5]) allow us algebraic exact or approximate 
analytical solution of ()3]) under some certain conditions. Let us substitute ([5]) into the ([I]) to 
obtain the function that is related to the wave function of the corresponding Hamiltonian: 

( in n \ / m n \ 

n=2 ' / \ n=2 ' / 

After all we have obtained an useful formalism of the Taylor expansion method. This form 
of the Taylor series can also be used to obtain series solution of the second order differential 
equations. In the solution of the eigenvalue problems, truncation of the the asymptotic 
expansion to a finite number of terms is useful. If the series optimally truncated at the 



smallest term then the asymptotic expansion of series is known as superasymptotic 13] , and 
it leads to the determination of eigenvalues with minimum error. 

Arrangement of the boundary conditions for different problems becomes very important 
because improper sets of boundary conditions may produce nonphysical results. When only 
odd or even power of x collected as coefficients of /(0) or f'(0) and vice verse, the series 



is truncated at n = m then an immediate practical consequence of these condition for 
<?m-2(0) = or p m _2(0) = 0. In this way, the series truncates at n = m and one of the 
parameter in the q m -2^) or Pm-2(0) belongs to the spectrum of the Schrodinger equation. 
Therefore eigenfunction of the equation becomes a polynomial of degree m. Otherwise 
the spectrum of the system can be obtained as follows: In a quantum mechanical system 
eigenfunction of the system is discrete. Therefore in order to terminate the eigenfunction 
f\x) we can concisely write that 

9m(0)/(0)+p m (0)/'(0) = 

g m -i(o)/(o) + Pm _ 1 (o)/ / (o) = o (7) 

eliminating /(0) and /'(0) we obtain 

gm(0)p ro _i(0) - Pm (0)q m -i(0) = (8) 

again one of the parameter in the equation related to the eigenvalues of the problem. 

We can state that the ATEM reproduces exact solutions to many exactly solvable differ- 
ential equations and these equations can be related to the Schrodinger equation. It will be 
shown in the following section ATEM also gives accurate results for non-solvable Schrodinger 
equations, such as the sextic oscillator, cubic oscillator, deformed Coulomb potential, etc. 
which are important in applications to many problems in physics. This asymptotic ap- 
proach opens the way to the treatment of Schrodinger type equation including large class 
of potentials of practical interest. 



III. SOLUTION OF THE SCHRODINGER EQUATION BY USING ATEM 

An analytical solution of the Schrodinger equation is of high importance in nonrelativistic 
quantum mechanics, because the wave function contains all necessary information for full 
description of a quantum system. In this section we take a new look at the solution of the 
Schrodinger equation by using the method of ATEM developed in the previous section. Let 
us consider the following eigenvalue problem (h = 2m = 1) 

- -J^ + V(x)i/>(x) = Ei{j(x) (9) 

where V(x) is the potential, ip(x) is wave function and E is the energy of the system. The 
equation has been solved exactly for a large number of potentials by employing various tech- 



niques. In general, it is difficult to determine the asymptotic behavior of fl9]) in the present 
form. Therefore it is worthwhile to transform (Q to an appropriate form by introducing 
the wave function ij)(x) = f(x) exp ^— J W(x)dxj . Thus, this change of wave function 
guaranties lim ip(x) = 0. We recast and we obtain the following equation 

x— >oo 

L(x) = -f\x) + 2W(x)f'(x) + (V(x) + W'(x) - W 2 (x) - E)f(x) = 0. (10) 

In this formalism of the equation coefficients in ((3]) can be expressed as: 

p (x) = 2W(x), q (x) = (V(x) + W'(x) - W 2 (x) - E). 

Using the relation given in §5§ one can easily compute p n (x) and q n (x) by a simple MATH- 
EMATICA program. Our task is now to illustrate the use of ATEM to obtain explicit 
analytical solution of the Schrodinger equation including various potentials. 

1. Anharmonic oscillator 

Solution of the Schrodinger equation including anharmonic potential has attracted a lot 
of attention, arising its considerable impact on the various branches of physics as well as 
biology and chemistry. The equation is described by the Hamiltonian 

H = -l- 2+ x 2 + gx\ (11) 

In practice anharmonic oscillator problem is always used to test accuracy and efficiency 
of the unperturbative methods. Let us introduce, the asymptotic solutions of anharmonic 
oscillator Hamiltonian when W(x) = x, then the wave function takes the form 

ifj = e 2 f(x), 

and (lT0~l)can be expressed as 

d 2 f df 

L(x) = ~^ + 2x£ + (gx* + 1 — E)f = 0. (12) 
Comparing the equations (j3J) and ( fl2l) we can deduce that 

Po(x) = 2x and qo(x) = (gx 4 + 1 — E). (13) 
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Here we take a new look at the solution of the ffl2l) by using the method of ATEM developed 
in the previous section. By applying flH]), the corresponding energy eigenvalues are calculated 
by the aid of a MATHEMATICA program. 

The term asymptotic means the function approaching to a given value as the iter- 
ation number tends to infinity. By the aid of MATHEMATICA program we calcu- 
late eigenvalues .Eand eigenfunction f(x) for g — 0.1 using number of iterations k = 
{20,30,40,50,60,70,80}. The eigenvalues are presented in Table I and and are compared 

n n 

with results computed by the AIM [1] and direct numerical integration method [14J by taking 
= 0.1. 



k 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


20 


1.06529019 


3.30632658 


5.74394288 


8.30136568 


10.82233628 


15.77123716 


30 


1.06528554 


3.30688248 


5.74807461 


8.35361636 


11.10902943 


13.84616319 


40 


1.06528550 


3.30687176 


5.74795553 


8.35266513 


11.09828874 


13.97716619 


50 


1.06528550 


3.30687202 


5.74795940 


8.35267765 


11.09860503 


13.96951294 


60 


1.06528550 


3.30687201 


5.74795926 


8.35267786 


11.09859535 


13.96995159 


70 


1.06528550 


3.30687201 


5.74795926 


8.35267782 


11.09859562 


13.96992450 


80 


1.06528550 


3.30687201 


5.74795926 


8.35267782 


11.09859562 


13.96992632 


E[l] 


1.065286 


3.306871 


5.747960 


8.352642 


11.09835 


13.96695 


E [14] 


1.065286 


3.306872 


5.747959 


8.352678 


11.09860 


13.96993 



TABLE I: Eigenvalues of the (llip for different iteration numbers k and g = 0.1. Last two rows 

n 

corresponds the comparison of eigenvalues computed by the AIM [1], direct numerical integration 
method Q . 

The function f(x) for n = 2 state is given in ffl4|) . 
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FIG. 1: The first four wavefunctions for anharmonic oscillator. 



k = 30; f(x) = 1 - 2.37404x 2 + 0.147996a; 4 + 0.0193758x 6 

-1.73022 x 10"V - 5.18734 x 10~ 5 s 10 + 8.68491 x 10"V 2 

k = 50; f(x) = 1 - 2.37398a; 2 + 0.14797a; 4 + 0.0193735a; 6 

-1.73037 x 10~ 3 a; 8 - 5.19242 x 10" 5 x 10 + 8.67726 x 10" 6 x 12 (14) 

k = 80; f(x) = 1 - 2.37398a; 2 + 0.14797a; 4 + 0.0193735a; 6 

-1.73037 x 10~ 3 a; 8 - 5.19243 x 10" 5 a; 10 + 8.67725 x 10" 6 a; 12 

As we mentioned before using ATEM we can obtain an analytical expression for the wave 
function of the Schrodinger equation. Substituting E into ((6]) we get the wave function of 
the Schrodinger equation for the corresponding eigenvalues. Analytical expressions signal 
that ATEM produce an efficient result for the eigenf unction. For the first four states the 
plot of the normalized wave functions are given in Figure 1. 
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2. Interacting electrons in a quantum dot 



In this section we present a procedure to solve the Schrodinger equation of two interacting 
electrons in a quantum dot in the presence of an external magnetic field by using ATEM. The 
problem has been discussed in various articles 



15 



17J . Here we just solve the mathematical 



part of the problem. Without further discussion the Schrodinger equation for a quantum 
dot containing two electrons in the presence of the magnetic field B perpendicular to the 
dot is given by 

H = Bi (Pl + eA{n))2 + + Ti^T] ^ 

i=i 1 ' ' 

Introducing relative and center of mass coordinates r — Ti — r\, R — \{f\ + r^) the 
Hamiltonian can be separated into two parts such that H = 2H r + \Hr, where 

p 2 1 e 2 1 

H r = — — + -m*u 2 r 2 + h -uJ c L r , (16a) 

2m* 2 2er 2 

H R = + \m*u 2 R 2 + \u c L R . (16b) 
2m* 2 2 

Equation (116bl) is the Hamiltonian of the harmonic oscillator, and it can be solved exactly. 
Let us turn our attention to the solution of the Hamiltonian H r . In the polar coordinate 
r = (r, a), if the eigenfunction 

cf) = r -h iea u{r) (17) 
is introduced, the Schrodinger equation H r (f> = E r <fi, can be expressed as 

- — ^ + —U 2 - t)4 + -mVr 2 + — + -UcL^ u(r) = E r u(r). (18) 
2m* dr 2 2m* y 4 J r 2 2 2er 2 J v ' v ' v ' 

From now on we restrict ourselves to the solution of Eq. (1181) . After changing the variable 



r 



— >■ / A_ _ r and substituting u(r) = r e+ 2e /(r),we obtain the following equation 



L( r ) = -rf"{r) + (ur 2 - {21 + l))/'(r) - (rE n + A)/(r) = 0, (19) 

where 

1 e 2 
E r = E n + (\£\+l)hu+-L r u ct x = ~^ (20) 

for simplicity we have chosen that fuo = = 1. In this case the functions po(r) and qo(r) 
are given by 

p (r) = ur - 2 ^ 1 and q (r) = -E n - ^. (21) 



Meanwhile we bring to mind that Hamiltonian ( TT5|) possesses a hidden symmetry. This 



implies that the Hamiltonian is quasi-exactly solvable [15|, Il6j . Fortunately, quasi exact 
solvability of the Hamiltonian gives us an opportunity to check accuracy of our result and 
to test our method. In order to obtain quasi exact solution of f[T§|) we set in: 

E n = ju, where j — 1,2, 3.... 

and then the problem is exactly solvable when the following relation is satisfied: 

A = {±vM2^+l)};j = l 
{o,±V2u;(4£ + 3)};j = 2 

jo, ±yj 10u(e + 1) ± ujV73 + 128£ + 64^ 2 | ; j = 3 



Note that A,o> and I belong to the spectrum of the Hamiltonian. Therefore an accuracy 
check for the ATEM can be made. We have tested ATEM and the result are given by 



f(r) = ItJ — — — r 
w 2£+ 1 



A 



{±^(47+3)} ;j = 2;E n = 2u; 



^ . /2o;(4£ + 3) 
/( r ) = 1= F\/ , 1 r + 



2^+1 2^+1 

Consequently, we demonstrated that our approach is able to reproduce exact results for 
the exactly solvable second order differential equations. Let us turn our attention to the 
complete solution of the ( TI5]) . We have again used 80 iterations during the solution of the 
equation and controlled the stability of the eigenvalues. The results are given in Table II. 

We go back ([6]) to obtain the wave function of the equation of (fl9l) for various values of 
E. Their plots are given in Figure 2. 



IV. CONCLUSION 



The basic features of our approach are to reformulate Taylor series expansion of a func- 
tion for obtaining both eigenvalues and eigenfunctions of the Schrodinger type equations. 
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k 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


n = 5 


20 


0.63844692 


0.81078941 


2.43413075 


2.80439505 


4.30282729 


4.81597264 


30 


0.64590375 


0.80463157 


2.45078087 


2.79146457 


4.32486753 


4.79696796 


40 


0.65007015 


0.80121630 


2.45998729 


2.78436605 


4.33697615 


4.78670663 


50 


0.65279850 


0.79898546 


2.46598524 


2.77974956 


4.34484756 


4.78008225 


60 


0.65475424 


0.79738773 


2.47027173 


2.77645083 


4.35046813 


4.77536829 


70 


0.65624076 


0.79617357 


2.47352334 


2.77394739 


4.35473042 


4.77180017 


80 


0.65741800 


0.79521187 


2.47609490 


2.77196620 


4.35810109 


4.76898147 



TABLE II: Eigenvalues E n of the (fT9|) for different iteration numbers k and A = \fuj = 1 and 
t = 1/2. 



FIG. 2: The wavefunctions of the two electron interacting in the harmonic oscillator potential field. 
The parameters A = yfui = 1. 
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Furthermore the technique given here has been applied to determine the eigenvalues and 
the eigenf unctions of the anharmonic oscillator and the Hamiltonian of two electrons in a 
quantum dot. We have shown that ATEM gives accurate results for eigenvalue problems. 



As a further wor 
mass Hamiltonians 



AT 
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M can be developed in various directions. Position dependent 
20| can be solved by extending the method given in this paper. In 
particular, Lie algebraic or bosonic Hamiltonians can be solved within the framework of the 
method given here. Before ending this work a remark is in order. This extension leads to the 
solution of various matrix Hamiltonians, Dirac equation and Klein-Gordon equation. Our 
present results manifest that ATEM leads to the solution of the Schrodinger type equations 
in different fields of physics. 
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